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FOREWORD 


This report describes a portion of the results obtained on 
NASA Grant NSG 3044. This work was done under subcontract to the 
University of Illinois, Urbana, with Prof. S.S. Wang as the Prin- 
cipal Investigator. The prime grantee was the Massachusetts Insti- 
tute of Technology, with Prof. F.J. McGarry as the Principal Inves- 
tigator and Dr. J.F. Mandell as a major participant. The NASA - 
LeRC Project Manager was Dr. C.C. Chamis. 

Efforts in this project are primarily directed towards the de- 
velopment of finite element analyses for the study of flaw growth and 
fracture of fiber composites. The analysis of such problems using 
three-dimensional analyses may be limited by the presence of very 
localized, high stress gradients as occur at free edges. The work 
described in this report is a theoretical investigation of such ef- 
fects resulting from thermal or hygroscopic loading. The results 
given here are for thermal effects, but the same treatment also ap- 
plies to hygroscopic effects, differing only by the expansion co- 
efficient used. Analogous results to those given here, but for hy- 
groscopic loading, may be found in the following papers: 

(1) S.S. Wang and I. Choi, AIAA paper 80-0713-CP, 21st Proceedings 
of AIAA/ASME/ASCE/AHS Structures, Structural Dynamics and Mate- 
rials Conf., Seattle, WA, 1980; also to appear in AIAA Journal, 
November, 1982. 

( 2 ) S.S. Wang and I . Choi , in Developments in Composite Materials 
and Structures, J.R. Vinson, ed., ASME , p. 315 (1980). 



ABSTRACT 


Thermal stresses near geometric boundaries of fiber-reinforced composite 
laminates play an important role in controlling complex failure modes and 
ultimate performance of composite materials in severe thermal environment. 

This paper presents an investigation of boundary-layer thermal stress singularity 
and distributions in composite laminates. Based on the theory of anisotropic 
thermoelasticity, a system of coupled governing partial differential equations 
is obtained. Edge boundary conditions and interface continuity conditions lead 
to a transcendental characteristic equation for determining thermal stress 
singularity at the laminate boundary. Complete thermal boundary-layer stress 
and displacement solutions are obtained by an eigenfunction expansion method 
in conjunction with a boundary collocation procedure. The thermoelasticity 
solution in the region away from the singular domain is found in excellent 
agreement with existing approximate numerical results. As the edge is approached, 
the singular terms control the near-field behavior of thermal boundary-layer 
stresses. Results are presented for the cases of various angle-ply graphite-epoxy 
laminates with [0/-0/-0/0] fiber orientations. Thermal boundary-layer thickness 
is defined by considering strain energy density distribution along ply interface. 
The thermal boundary-layer thickness is shown to depend on the anisotropy of 
individual lamina, ply thermomechanical properties, and relative thickness of 
adjacent layers. 
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1. INTRODUCTION 

The response of a composite laminate near its geometric boundaries sub- 
jected to severe thermal and other environmental loading has attracted much 
attention recently, since failure of composite materials is frequently initiated 
at the boundaries. The thermal stress field in the vicinity of the laminate 
boundaries, which may be primarily responsible for strength degradation and 
failure of composites, has been investigated by several researchers using 
different approximate methods [1-4] . Approximate solutions for the thermal 
boundary-layer problem have revealed several unusual features. Interlaminar 
stresses near a traction-free boundary of a composite laminate subjected to a 
uniform temperature change have been found to be very high and inherently 
three-dimensional. It has also been reported that the high thermal stress is 
confined within a localized region of several laminar thicknesses from the 
edge and that in the boundary-layer region they cannot be assessed accurately 
by the classical lamination theory [5,6]. The behavior of this highly stressed 
boundary-layer region is of great importance in controlling complex failure 
modes and ultimate performance of the composites. Understanding the fundamental 
nature of boundary-layer thermal stresses is essential to the failure analysis, 
design, and processing optimization of composite materials. 

While all previous numerical approximate solutions indicated an unbounded 
trend of boundary-layer thermal stresses and postulated possible existence of 
a stress singularity at the edge, the search for the order or strength of the 
stress singularity has been unsuccessful. No information concerning the exact 
order of the boundary-layer stress singularity has been reported yet, to the 
authors^ knowledge. The apparent difficulty may result from the complicated 
nature of the problem such as the thermomechanical anisotropy of each individual 
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fiber-reinforced lamina, the geometric discontinuity, and the abrupt change of 
material properties through the laminate thickness direction. Since the 
thermal boundary-layer effect is localized in nature, it is necessary to 
determine the exact order of the laminate edge stress singularity so that 
complex thermal response in the vicinity of laminate boundaries can be studied 
more accurately. This paper presents a rigorous theoretical study of thermal 
boundary-layer stress singularity and distributions in composite laminates 
subjected to uniform thermal loading. 
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2. FOEMULATION 

2.1 Basic Equations 

Consider a general composite laminate composed of fiber-reinforced 
plies subjected to mechanical and thermal loading, as shown in Fig. 1. Denote 
the constitutive equation of each individual ply by the Duhamel-Neumann form 
of generalized Hooke’s law in contracted notation as 

G.=S,.a.+a, AT, (i,j=l,2,3,....6), (1) 

X XJ J 1 ’ J J 


where the repeated subscript indicates summation; is the compliance tensor; 
a^, the thermal expansion coefficient, and AT, the change of temperature. The 
engineering strains, e^, in Eq. 1 are defined in a Cartesian coordinate system 
by 




Go = 


z 

z 


3w 

3z’ 


o , 9v 

e, = 2e = 

4 yz 9y 9z 


9w . 9u 

G- = 2e = — + — , 
5 xz 9x 9z 


- 9u 9v 

G^ = 2g = 

6 xy 9y 9x 


( 2 ) 


where u, v and w are displacement components. The stresses, are 
defined in an analogous manner. 

The composite laminate considered here has a finite width and is sub- 
jected to a uniform axial extension, e, along the z-axis and a uniform temper- 
ature change, AT. The composite is assumed to be sufficiently long that, in 
the region far from the end, the end effect is neglected by virtue of Saint 
Venant’s principle. Consequently, stresses in the laminate are independent 
of the z coordinate. The special case where e vanishes identically corresponds 
to the well known generalized plane deformation [7]. Under these assumptions, 
equations of equilibrium without body force read 


9a 

X 

9x 


+ 


9t 


9y 


0 




(3a) 
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Integrating Eq. 


u = 


V = 


w = 


dT 3a 

4. _Z ^ 

3x 3y 


0, 


9t 3t 

-^+-^= 0 . 

3x 3y 


1 with the aid of Eq. 2, one can 
z2 9D 

2 3x 5j j 5 3x 

2 

- ^ 1^ + (S, a. + a AT - -^) 
2 9y 4j j 4 3y 


D z + W^(x,y), 


obtain u, v and w as 
z + U^(x,y) , 

z + V^(x,y) , 


(3b) 


(3c) 


(4a) 

(4b) 

(4c) 


where 

D = o, + AT, (5) 

and U , V and W are arbitrary functions of x and y only. 

Following the procedure in [7], it can be shown after some mathematical 
manipulation that general expressions for the displacements and the stress 
component have the following forms; 


U = - ^ z 2 - yz + U(x,y) + (d^ z - o)^ y + u^, (6a) 

V = - Y ^2 ^33 xz + V(x,y) + x - z + v^, (6b) 

w = (A^ X + A^ y + ^2)8^2 z + w(x,y) + 03^^ y - o)2 x + w^, (6c) 

X + y + A3 -(83^ Cj + «3 AT)/S33, (j = 1,2, 4 , 5 , 6). (6d) 


The unknown functions, U, V and W, depend on x and y only, and can be shown 
easily to obey the following relationships: 


3x 



a. + 
3 


SiaCAi X + y + A^) + AT 


(7a) 
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3V ~ 

^ = S 2 J Oj + X + A 2 y + A 3 ) + 02 AT, (7b) 

3W ~ 

^ = Sjj 0 ^ + S 33 <Aj^ X + A 2 y + A 3 ) + A^ y + o^ AT, (7c) 

3W 

87 = % + ^A3^\ X + A 2 y + A 3 ) - A^ X + o^ AT, (7d) 

37 37 - ®6J + ®63<*1 * + *2 y * -6 


where 


(j = 1 , 2 , 4, 5 , 6). 


^ij ^ij ■ ^13 ^j3'^^33’ 


( 8 a) 


“i = “i - “3 ^ 13/^33 


(i,j = 1,2, 4, 5, 6). (8b) 


It is obvious that the constants, u^, v^, and = 1>2,3), in Eqs. 6a-6d 

characterize the rigid body translation and rotation of the solid. and A 2 

represent the bending of the laminate in the x-z and y-z planes, respectively. 
A^ characterizes the uniform axial extension of the composite laminate, and A^, 
the relative angle of rotation about the z-axis. 

2.2 Governing^ Partial Differential Equations^ and Boundary and End Conditions 
Introducing Lekhnitskii^s stress functions, F and such that 


a 

X 


~3y^* 


8^ F _ 8^ F 

^y 3x^ ’ ^xy 9x 3y* 


= 11 
^xz 9y* 


3y 

^yz 3x* 


(9a-e) 


one can show that the equations of equilibrium are satisfied identically. 
Eliminating U and V from Eqs. 7a, 7b and 7e, and W from Eqs. 7c and 7d, one 
obtains the following system of governing partial differential equations: 
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fLj F + Lj f - - 2 A^ + S34 - Aj S35 - AT, 


F + 4' = - H 2 AT, 


(10a) 

(10b) 


where L 2 , L^, L^, and H 2 are linear differential operators defined as 


= S 


- 2S 




'2 44 3x2 ^^45 3x 3y ^ ^55 3y2» 


(11a) 


"3 ^^24 ^^25 ^46^3x2 3y “ ^^14 ''' ^56^3x 3y2 ^15 3y3‘ 


L. = - S. 


S ^22 3x*» ^^26 3x3 3y ^^^12 ®66^3x2 3y2 


- 2S 


a4 ~ 3 ^ 

^ + S, ^ 


16 3x 3y3 “'ll 3y4» 


and 


»i = -“4lr+“5^’ 


**2 “ “2 “6 “1 W ' 


(11c) 

(lid) 

(lie) 


Now consider boundary conditions on the lateral surfaces and at the ends 
of the laminate strip. Assuming that the lateral surfaces are free from sur- 
face tractions, one may have the following conditions: 


and 


a n + T n 
XX xy y 


T n + 0 n 
xy X y y 


T n + T n 
xz X yz y 


= 0, 
= 0, 

= 0, 


(12a) 

(12b) 

(12c) 


where n and n denote the directional cosines of the outward unit normal on 
X y 

as shown in Fig. 2. The conditions at the ends of the composite laminate 
may be obtained from the statically equivalent loading as 
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n 

If 

] 

II 

II 

] 

II 


X dx dy = 0 
xz. 


X dx dy = 0, 
yz ’ 


a dx dy = P , 
z z 


y dx dy = M^, 


a X dx dy = M , 
z y 


B 

I (t X - t y) dx dy = 

JJ yz xz ■' t’ 

B 


(13a) 


(13b) 


(13c) 


(13d) 

(13e) 

(13f) 


at z = ±L, where P , M , M , M are the applied axial force, bending moments, 
z X y X 

and twisting moment , respectively. 


2.3 Interface Continuity Conditions 

Consider a portion of the laminate composed of two different fiber- 
reinforced laminae, as shown in Fig. 2. Stress functions of Eqs. 9a and 9b 
must be defined for each ply. Assuming that the plies are perfectly bonded 
along the interface one can immediately establish the continuity con- 

ditions of the stress and displacement along the interface between the kth and 
(k+l)th plies as the following: 


^(k) ^(k) . ^(k) ^(k) ^ ^(k+1) ^(k+1) ^(k+1) ^(k+1) 

an+xn=-a n -x n , 

xxxyy XX xyy 

(k) (k) (k) (k) (k+1) (k+1) ^(k+1) (k+1) 

xy X y y xy x y y 

(k) (k) ^ (k) (k) (k+1) (k+1) (k+1) (k+1) 

xz X yz y xz X yz y 


(14a) 

(14b) 


(14c) 
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u(k) ^ ^(k+1)^ 

v(k) = 


(15a) 

(15b) 

(15c) 


By using divergence theorem, it can be readily seen that Eqs. 13a and 13b are 
satisfied identically by virtue of 3c, 12c and 14c. Substituting the displace- 
ments, Eqs. 6a-c, into the above continuity conditions, Eqs. 15a-c, one finds 
that the unknown constants for two adjacent layers in Eq. 6 are related by 


(k) _ , (k+1) ^ (k) _ ^ (k+1) 

o o * o o 


(k) (k+1) 

w = w 
o o 


(16a-c) 




(k) _ (k+1) 

“2 ■ “2 


(16d-e) 




and 


,(k) „(k) ^ ,(k+i) „(k+i) 
i 33 i 33 


(i = 1,2,3) 


^(k+1)^ (16f-h) 


(16i) 


^(k) ^ ^(k+1) 
4 4 


(16 j) 
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3. METHOD OF SOLUTION 

The complex governing differential equations and boundary conditions 
formulated in the previous section may be simplified considerably as 
appropriate loading conditions and geometric symmetry are taken into con- . 
sideration. For simplicity and without loss of generality, the following 
geometric and environmental conditions are introduced for the present thermal 
boundary-layer stress problem: 

(a) the temperature change AT is constant and uniformly distributed 
throughout the laminate, and there is no external mechanical loading applied; 

(b) the interface is a straight line and meets the traction-free edge of 
the composite laminate by a right angle. 

Under these assumptions the governing differential Eqs.lOa and 10b may 
be simplified as 

|L 3 F + T = - 2A^ + - A 2 S 33 , (17a) 

< 

F + L 3 <P = 0, (17b) 

with boundary and end conditions, 

CT=x=x= 0 , onx= 0,(18a-c) 

X xy xz 


II dx dy = || x dx dy = || 0 ^ y dx dy = 0, on B, 


(19a) 


B 




on B. (19b) 


B 


The interface continuity conditions of tractions and displacements may be ex 
pressed as the following: 
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^(k) ^ ^(k+1) 


^(k) _ ,(k+l) 

T = T 


xy 


xy 


T<« = T»«>, (20a) 
yz yz 


y(k) ^ u(k+l) 


X = V 


(k+l)^(k+l) 


^(k) ^ y(k+l) (20b) 


along y = 0, as the original of the coordinate system is moved to the free edge, 
3.1 ■ Eomogeneous Solution 

The simplified governing differential equations, Eqs. 17a and 17b, are 
coupled, linear partial differential equations with constant coefficients re- 
lated to the anisotropic thermomechanical elastic constants of each composite 
lamina. The general solution consists of a homogeneous solution and a par- 
ticular solution depending on geometric, loading, and boundary conditions of 
the problem under consideration, Lehknitskii has shown [7] that the homo- 
geneous solution of the above mentioned system of governing partial differential 
equations has the general form 


F(x,y) = \ F (x + 


(21a) 


k=l 


’l'(x,y) = I Fj^(x + y). 


k=l 


(21b) 


where the prime (*) in Eq. 21b denotes differentiation of the function 
F^(x + y) with respect to its argument, and the coefficients \x^ are the 
roots of the following algebraic characteristic equation: 


and 


with 


£^(m) ^2(y) - ^3(y) = 0. 




(22a) 

(22b) 


^55 ^^45 ^ ^ ^ 44 ’ 


(22c) 
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^3^^^ = S^5 ~^^14 '^^^25 ■*■ ^46^ “ ^24’ 

^4(m) = y*t - 2S^g y3 '^^2^12 ^66^ " ^^26 ^ ''' ^22* 

It has been shown [7] that Eq. 22a cannot have a real root. Thus the roots 
appear as complex conjugates, and are analytic functions of the complex 
variables = x + y. Substituting the expressions of F(x,y) and 't'(x,y) 
in Eqs. 21a and 21b into Eqs. 9a-e, the homogeneous components of the stress 
may be expressed in terms of 


X 

- j, < W- 

k=l 

(23a) 

y 

6 

= I 

k=l 

(23b) 

(h) 

"^yz 

^ ■ X 

k=l 

(23c) 

^(h) 

xz 

■ J, \ \ 

k=l 

(23d) 

,(h) 

xy 

- - X -k 

k=l 

(23e) 


The expressions for displacement components may be obtained directly from 
Eqs. 7, 10 and 23 with omission of terms which are to be included in 
the particular solution. Hence, one has 


- I »k w- 


k=l 




k=l 


(24a) 


(24b) 
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and 



(24c) 

where 

Pk 

■ ®U + ®12 - ®14 \ + ®15 \ \ - ®16 “k> 

(24d) 


\ 

®12 “k '*' ®22^\ ■ ®24 \^‘‘k '*' ^25 "k " ^26’ 

(24e) 


\ 

' ®14 \ ■" ®24''“k - ^44 "k^“k ®45 "k - =46- 

(24f) 

We 

now choose the form of 




S Z^'^2/[(6+1)(5+2)], 

(25) 

where C, 
k 

and 5 

are arbitrary complex constants to be determined later. 

Sub- 


stituting Eq. 25 into Eqs. 23 and 24 gives 


and 


X 


a^^> = 

y 



^(h) 

xz 



, ^k+3 ^k 

k=l 

J ‘‘ic 4 + \+3 

k=l 

- J t<=k \4 + "k+3 

k=l 

3 . , 

J,“=k \ “k \ + \+3 \ \ 

k=l 

“ ^ ^k \ ^k+3 

k=l 


- . I ""k »k C'" + V3 »k " 


,5+1 


t;5+1. 


k=l 


(26a) 


(26b) 


(26c) 


(26d) 


(26e) 


(27a) 
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= X “"k "k 4*’’ * Vs S - (27») 

k=l 

k=l 

where the overbar denotes the complex conjugate of the associated variable. 

The homogeneous solutions for the stress and displacement shown in 
Eqs. 26 and 27 are required to satisfy homogeneous boundary conditions 
and interface continuity conditions, i.e., Eqs. 18 and 20. This leads to a 

standard eigenvalue problem, which results in a complicated transcendental 
equation for determining the value of 6. By taking the complex conjugate 
of Eqs. 26 and 27, It is readily seen that 6^ always appear as pairs of 
complex conjugates. Thus the stresses and displacements can be made real by 
superposing the conjugate solutions. It is noted that a similar method of 
extracting the real function for an isotropic material case was discussed by 
Theocarls [8] who assumed the stress function including both 6 and 6 in the 
formulation. It is also noted that, due to positive definiteness of 
strain energy of an elastic body, the value of 6^ bounded by -1 < Re [6^] < 0 
characterizes the order of singularity of the boundary-layer stress field in 
the composite laminate. It has been found in a related paper [9] that there 
is only one real 6^ which meets this condition in the singular boundary-layer 
stress problem. Detailed information on the edge-stress singularity has been 
given in Reference 9. 

After the eigenvalues, 6^, are determined, the stress and displacement 
may be expressed in terms of the resulting eigenfunctions, f^^ and 

^ V = 1, 2,4,5, 6), (28a) 

m 
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,(h) _ 


I d g.„(x,y; 5), 


m 


(j = 1,2,3), 


(28b) 


where the unknowns d are real, and f. and g, denote known eigenfunctions 

m * im ^ 

corresponding to the m-th eigenvalue 6^. 

3.2 Partiaular Solution 

The particular solution for the problem may be sought in the form of 
polynomials as 

(x,y) = a^ +a^ y i-a^ xy^ + a^ y^ + a^ x^ + a^ xy + a^ y^ , (29a) 


^^^^(x,y) = ag x^ + ag xy + a^^ y^ +a^^^ x + a^^ y» 


(29b) 


where a^ are arbitrary constants to be determined. It is seen that Eq. 17b 
is satisfied identically and that Eq. 17a provides the following relation for 


each ply: 




+ + 2S^^ ag - 2S^^ a^ + 2S^^ a^^ 


- -2^4 + - Ag Sgg , 


Substituting Eqs. 29a and 29b into Eqs. 7 and 9 gives 


= 2a„ X + 6a, y + 2a_, 
X 3 4 *' 7 


(31a) 


= 6a^ X + 


= 6a- X + 2a« y + 2a_, 
y 1 2 3 


(31b) 


(p) 

^yz = - 2ag X - ag y - a^^^. 


(31c) 


^xz^ = ^9 ^ + ^^10 y ^12’ 


(31d) 


= - 2a„ X - 2a y - a., 
xy 2 3 6 


(31e) 
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3U 


(p) 


ax + + 


av 


(p) 


ay = ^21 ^ + ^22 y + ^23’ 


aw 


(p) 


ay = ^^41 - ^ ^42 y + ^43’ 


aw 


(p) 


ax ^51 ^ ^^2 y ■*■ ^53’ 


where 


ay ax ^61 ^ ^62 ^ ^63’ 


^jl ^3 ^^j2 ^1 “ ^^j4 ^8 ^j5 ^9 ~ ^2 '*' ^j3 ^1’ 

^j2 = ^ + 2 Sj 2 a^ - ag + 25^^ - 2 S.g a 3 + 8^3 A^, 

^j3 = 2S.3L + 28.2 S " ^j4 ^1 ^j5 ^12 " ^j6 ^6 + ^j3 ^3 “j 

Integrating Eqs. 32a-e with the aid of Eq. 30, one can obtain 


U 


- E,i x 2 + E32 xy + X + i(E^, - E,,) y^ + ^ E^, y. 


2 11 


13 2' 62 21 ' 


2 63 


V^P^= xy + i E,, y2 + E,, y + ^(E^^ - E,,) x^ + ^ E^, x. 


23 


2^ 61 12^ 


2 63 


W^P>= J E33 x2 + (E32 + A^) xy + E33 X + f- E^2 + ^^3 y- 


Thus the particular solution for the displacement can be written as 

u^P^= - 833 z2 - A^ yz + U^P^(x,y) + 0)3 z - (1)3 y + u^, 

v^P^ ^ ^ ^2 ^33 + A^ xz + V^P\x,y) + 0)3 x - a)^ z + v^. 


/P)= 


(P) , 


(A3 X + A2 y + A3) 833 z + W ^ (x,y) + y - 0)3 x + 


(32a) 

(32b) 

(32c) 

(32d) 

(32e) 

(33a) 

(33b) 

(33c) 

(34a) 

(34b) 

(34c) 

(35a) 

(35b) 


(35c) 
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Equations 31 and 35 are required to satisfy the boundary conditions, Eqs. 18a-c, 
and the interface continuity conditions, Eqs. 20a-b of the current problem. 

This leads to the establishment of the following relations 



2 hz 



p(k+l) 

®63 




(k+1) (k) _ p(k+l) 

51 ’ ^3 " ^3 • 


(36j-l) 


By examining Eqs. 16, 30 and 36, it is observed that there are 44 unknowns (in- 
cluding the unknown rigid-body translations and rotations) related by 34 linear 

algebraic equations. Solving these equations, there remain ten unknowns (for 
(k) (k) 

example, ), which may be determined by Eqs. 19a-b and 

the boundary conditions other than those on the traction-free edges. 


3 . 3 Complete Solution 

Now the complete solution for the thermal boundary-layer stress problem 
can be written as 


(h) 

a± = ^ 

1 

(1 = 1,2, 4,5, 6) , 

(37a) 

(h) 

u. = u. 

3 J 


(j = 1,2,3) , 

(37b) 


where given by Eqs. 28a-b and Eqs. 31, 35, 

respectively. The solution for the stress, Eq. 37a, satisfies identically 


the boundary conditions, Eqs. 12a-c, of the free-edge surface 9B . Along 

F 

> residual stresses appear due to the particular 
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solution being posed in the previous section. The residual stress can be 
counterbalanced by the homogeneous solution; thus, on one obtains 


Xh) 

n 

+ 

,(h) 

n = 

-a<P> 

n 

-x(P> 

n , 

(38a) 

X 

X 


xy 

y 

X 

X 

xy 

y 

(h) 

n 

+ 


n = 


n 


n , 

(38b) 

xy 

X 


y 

y 

xy 

X 

y 

y 

(h) 

n 

+ 

,(h) 

n = 

-x^P> 

n 


n . 

(38c) 

xz 

X 


yz 

y 

xz 

X 

yz 

y 


Without orthogonality among eigenfunctions, Eqs. 38a-c may be satisfied 
numerically in a least square sense through a boundary collocation method 
by truncating the eigenfunction series. It is noted that the particular 
solution is coupled with the homogeneous one through the end conditions, 

Eqs. 19a-b. Thus by matching Eqs. 38a-c with the aid of Eqs. 19a-b, one can 
determine all the unknowns explicitly. To this end, the expressions for 
may be obtained from Eqs. 6 d and 37a as 

= - S 3 J /S 33 (j = 1 , 2 , 4, 5, 6 ), ( 39 a) 

X + A 2 y + A 3 - ( 83 ^ + 03 AT)/S 33 . (39b) 

Complete expressions for the displacement field can be written explicitly in 


a similar manner. 
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4 . NUMERICAL EXAMPLES AND DISCUSSION 


For Illustration, symmetric angle-ply composite laminates with [±6] ori- 

s 

entations are studied. Thermal loading in the form of a uniformly distributed 
temperature change is considered. The particular laminate configuration and 
fiber orientation are chosen because approximate numerical solutions for this 
problem are available in the literature for comparison. Each individual lamina 
is assumed to be high-modulus graphite-epoxy with the following elastic constants: 
Ell = E 22 = 2.1 X 1q6 (psl), E^2 = 20 X 1q6 (psi) , = ^13 = G23 = 0.85 x 10 ^ 

(psi), v ^2 = ^13 ’^93 = 0.21, = U 2 = 16 X 10~® (/°F), = 0.2 x 10 ^ (/“F), 

where the subscripts, 1, 2 and 3, refer to transverse, thickness, and longitu- 
dinal directions of the individual ply, respectively. The and are 
evaluated by using these material constants, and the following relationships 


can be established readily: 

= 0 (3=1, 2, 3, 5, 6); = 0 (i=l,2,3,5; 3=1,2); 

(i,j i 3); (j=4,5,6); (i,j ^ 4 and 1 j); 

0=1.2); 0 (6=1,2); 5f>=Sf>; if>=-jf>. (4( 


4.1 Symmetry Conditions and Further SimpZications 

The geometric and lamination symmetry conditions lead to the following 


relationships: 

U(x,y) = U(x,-y), V(x,y) = -V(x,-y), W(x,y) = W(x,-y), (41a-c) 

U(x,y) =-U(-x,y), V(x,y) = V(-x,y), W(x,y) =-W(-x,y). (42a-c) 

Equations 41a-c and 42a-c may be written in equivalent forms as 

U ,,(0,y) = V ^(0,y) = W ^ C0,y) = 0, (43) 

>y ,y 

U (x,0) = V (x,0) = W (x,0) = 0. (44) 

,y^ * ^ ,x^ ’ ^ ,y^ ’ ^ 


The relations provided by Eqs. 7 and 40 and the symmetry conditions given in 
Eqs. 41 and 42 suggest that 
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= 0 


(i=1.2,4), 


and that a^, a^, and are symmetric with respect to the x and y axes 
and T and x are antisymmetric with respect to the x and y axes. Thus, 

(b) 

only ' remains to be determined and only one-quarter of the laminate needs 
to be considered. 

By using Eqs. 40a and 40b, one can easily show that 


^5 - ^5 


' (1) 

~(1) 


+ ^ 

;(i) *3 


®52 

^52 


- m ( 2 ) 

“3 " “3 ’ 


A3 A^ 


and all other unknowns vanish. Thus, the particular solution takes the 
following forms: 


a'P'S - 0, 

X 


53 (1) 5 

~(1) ^3 “ -(1) 

52 52 


^(P)6 . ^(P)6 . 0, 
xz xy 


(49a) 


(49b) 


(49c) 


in both layers (3 = 1,2). The displacements may be shown to have the expressions, 

T „(i) ~(i)l f ~(i) 1 

u(p)3= Af>+ AT X, (50a) 

li ^52 J i ^2 J J 
T f 1 

„(P)3 ^ e(l) _ ,5_3_ ^2 2_ ,(1) ^ AT „ ..OKA 
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and 


= 0. (50c) 

Inserting Eqs . 49a-c into Eq, 39b, one can obtain the stress component 


as 

z 


,(P)3 _ 


1 + 




^32 

^3 

~( 1 ) 


^52 

^33 

J 


A 3 + 


^32 “5 


AT 

~(1) 

“3 


^52 

J 

^33 

may be established: 


(51) 




s(i) 

s^3)l 


dx dy = 4 
z 

1 + 

^53 


^ ~(1) 


0 i 


^52 

33 


h(h^ + h^) A^^^ 


B 


+ 4 


^32 “5 (1) 

~(1) " “3 

^52 


b(hi + h^) ^ 

33 


(52) 


The contribution from the homogeneous solution may be obtained from Eqs . 28 
and 39a by 


II ^ ■>*'<!' = I 1'3„' 

B ^ 


where 


3m 




(53) 


(54) 


Substituting Eqs. 52 and 53 into Eqs. 19 and solving the resulting equation 


for A„ in terms of d give 
3 m 


= e + I d e 

< n ^ 1 

m 


3 '•o ■ ^ m 3m’ 


(55) 


where ® 3 jjj^re known constants obtained by integrating associated eigenfunctions. 


and e^, from the particular solution. 
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The boimdary conditions on (8B - 3B ) suggest, by using Eqs . 28, 38, 43, 

r 

44, 49 and 50, the following relationships for determining the constants d : 

in 


and 


Id 


m 


m 


I d 
I d. 

m 

I d 


m 


m 


m 


m 


I d 

m 


m 


I d 

m 

I d 

m 

I d 


m 


m 


m 


m 


,( 2 ) 

2m 


(x,h ; 6 ) 
X m 

(x,hi; 6^) 

(x,hi; 

(x,-h2,6m) 

(x,-h.;6 ) 
z in 

(b,y; 6. ) 

m 

(b,y; SJ 


!53_^(i) 

:( 1 ) ^3 ^ :( 1 ) 


52 


AT, 


52 


0 , 


0 , 

0 , 

0 , 

0 , 


0 , 

0, (3 = 1,2) 


(56a) 

(56b) 

(56c) 

(56d) 

(56e) 

(56f) 

(56g) 

(56h) 


J <561) 

m 

where h. , denote differentiated forms of g , . according to Eqs. 43 and 44 with 
ij ij 

the origin of ' the Cartesian coordinates being transferred to the left 

free edge (Fig. 2). The boundary conditions shown in Eqs. 56a-i are matched 

by a boundary collocation method, using the eigenfunctions in Eq. 28. The 

constants d are then evaluated by satisfying these conditions at a given 
m 

number of selected collocation stations. Accuracy and convergence of 
solutions and effects of collocation points along the boundaries are 
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reported elsewhere [9]. For the numerical results reported in this paper, 

63 terms in the eigenfunction series and 74 collocation stations along the 
boundaries were used to ensure the accuracy and convergence of the solutions 
[9]. 

In what follows, laminate thermoelasticity solutions determined from 
the current eigenfunction expansion method are presented first and compared 
with existing approximate numerical solutions available in the literature. 
Detailed results characterizing the thermal boundary-layer field in the 
symmetric angle-ply composites with various lamination variables are given 
also. 

4.2 ThermaZ BouTidavy -Layer Stress Distrihution 

Distributions of the in-plane and interlaminar thermal stresses, a^, 

T » a and t , along the ply interface of a [45®/-45°/-45®/45®] graphite- 
xz y yz 

epoxy laminate are shown in a semi-logarithmic plot in Fig. 3. Dotted lines 
in the figure represent the result given in [3] by using constant-strain tri- 
angular elements (CSX) in a finite-element procedure. Solid lines denote the 
present laminate thermoelasticity solution with free-edge stress singularity 
being included. The two solutions are in good agreement in the region away 
from the laminate boundary. The in-plane stresses, and '^^^9 in the region 
away from the edge are found to be relatively constant and to recover to more 
or less what classical lamination theory (CLT) predicts. (Based on the 
classical lamination theory, the only stress induced by a unit change of 
temperature in the [±45]^ graphite-epoxy composite is a constant in-plane 

shear stress = 28.8 psi/°F.) As the edge is approached, the difference 
xz 

becomes gradually appreciable due to the presence of the stress singularity, 
which was not included in the previous approximate numerical solution. As 
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will be shown later, within the boundary-layer region the stress field is 
completely governed by the singular terms in the stress solution, and the 
interlaminar stress increases very rapidly. In fact, the interlaminar stress 
is much higher than the in-plane stress component as the edge of the 
laminate is inf initesmally approached. Thus thermally-induced deformation 
and failure may be dominated by the interlaminar stresses. 


4.3 Stress Singularity in Thermal Boundary -Layer Field 

Since the stress and strain fields within the boundary-layer region are 
governed by the singular terms in the present laminate thermoelasticity solu- 
tion, theoretically, the thermally induced stresses determined from the solution 
are unbounded at the intersection of the ply interface and the laminate edge. 
Thus the near-field stress may be expressed in a general form as 


a 


i 



+ 0 (higher-order, non-singular terms) 


(i=l,2,3,...,6) , (57) 

where Z and Z have their origin at the intersection of the ply interface and 
the edge of the laminate (Fig. 2); 6^ is the order of the thermal boundary- 
layer stress singularity which is the smallest eigenvalue satisfying 


0 > Re[6^] > -1 (58) 

among all the 6^ determined from the characteristic equation [9] in solving 

for the homogeneous solution of the governing partial differential equations. 

The order of the boundary-layer stress singularity is noted [9] to depend 

only upon lamina constitutive properties and fiber orientations of the adjacent 

plies. Numerical results of the first four nonzero eigenvalues 6 for sym- 
^ m 

metric angle-ply [0/-0/-0/9] graphite-epoxy laminates are given in Table 1. 
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Note that corresponds to the order of boundary-layer stress singularity 
and that zero and intergers, n, are always eigenvalues for the problem. In 
this composite system, 6^ has a value which is rather weak as compared with 
other typical singular stress problems such as an elastic crack problem. 

Among various 0 studied, it appears that the [±51®]^ graphite-epoxy laminate 
has the strongest stress singularity [9]. 

4.4 Thermal Boundary -Layer Stress Intensity Factors 

For a composite laminate with given fiber orientations, the coefficients 
of the singular terms in Eq. 57 characterize amplitudes of the thermal stress 
and strain in the boundary-layer region. Because the boundary-layer stresses 
are most crucial along the ply interface, i.e., the x-axis, and become singular 
at the interface/edge intersection, it is possible to define the amplitudes 
of the singular thermal boundary-layer stresses by 

—6 1 

K. = lim X (i=l,2,3, . . . ,6) . (59) 

^ xr>0 ^ 

The are dependent upon geometric variables of the composite (e.g., ply 
thickness, number of layers), lamination parameters (e.g., fiber orientation, 
stacking sequence), and mechanical and thermal loading conditions. The 
fundamental structure of the thermal boundary-layer stress solution shown 
in Eqs. 57 and 59 resembles that of an elastic crack problem (except that 6^ 
has a value of -0.5 in the crack-tip stress field). Also, the nature of 
is similar to the so-called crack-tip stress intensity factors in linear 
elastic fracture mechanics. Thus, in this context, it may be appropriate to 
denote as “thermal boundary-layer stress intensity factors" or "thermal 
free-edge stress intensity factors" for the composite laminate. Values of 
for the [0/-0/-0/0] graphite-epoxy composite with all laminae being of 
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equal thicknesses under unit thermal loading are determined in Table 2, in 

- 6 . 

which carry the unit of [psi-in.] 

The associated with the interlaminar stresses are found, in general, 
significantly larger than those associated with the in-plane stress components. 
The dominance of the interlaminar thermal stresses a and t in the boundary- 

y yz 

layer region shown in Fig. 3 is clearly illustrated by the high values of K 2 
and K^, which are, in fact, one or two orders of magnitude higher than the 
associated with a , a and x for all 0 studied. The high negative value of 

Z XZ 

indicates that a large compressive interlaminar normal stress is 

developed near the edge. Note that is found to vanish for all angle-ply 

[±9] composite laminates, due to the symmetry of ply orientations and traction- 
s 

free edge conditions. 

4.5 Through-Thickness Distribution of Thermal Boundary-Layer Stresses 

The unique features of the thermal boundary-layer effect are further 

illustrated by through-the-thickness distributions of in-plane and interlaminar 

thermal stresses near the laminate boundary. The in-plane thermal stress 

in the thickness direction at different distances away from the edge are shown 

in Fig. 4. The current solution is in agreement with previous results [3] 

that is compressive in the most part of the section near the laminate 

boundary except for the region closest to the ply interface where boundary- 

layer stress singularity dominates. Through-the-thickness distributions of 

the most dominant thermal interlaminar stress t are given in Fig. 5. The 

yz 

gradient of in the y-direction increases rapidly as the laminate edge is 
approached. Again the laminate elasticity solution in the boundary-layer 
region differs from the approximate solution [3] near the interface but in 
in good agreement in the far field. The next two figures (Figs. 6 and 7) 
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provide information on through-the-thickness distributions of thermal inter- 
laminar shear and normal stresses, and in the [45°/-45°/-45®/45°] 
graphite-epoxy laminate. The is found to change significantly with the 
thickness coordinate, alters its sign across the interface, and reaches a 
maximum value at a small distance y away from the interface. Within the 
boundary-layer region, t becomes very small along the interface and vanishes 
identically at the laminate boundary as envidenced by = 0. The distribution 
of at several distances near the boundary of the laminate is shown in 
Fig. 7. At any given x/b, is generally very small at a distance y away 
from the interface, and reaches a higher level as the interface y = h is 
approached. As one moves towards the edge, the Interlaminar normal stress 
becomes very significant in compression. As x 0 and y h the stress 
solution becomes unbounded due to the stress singularity at the intersection 
point. 


4.6 Thermal Boundary-Layer ^idth 

The rapid increase of thermal stress has been observed to be restricted 
to within a very localized region near the edge of the laminate — the so-called 
"thermal boundary-layer width." The thermal stresses developed in the 
boundary-layer region are inherently three-dimensional and cannot be determined 
by classical lamination theory. The singular nature and the extent of pertur- 
bation of the thermal boundary stress are considered to be of vital importance 
in controlling initiation of interlaminar fracture (or delamination) and 
strength degradation. The extent of perturbation of the laminate thermal 
stress field can be characterized by the boundary-layer width (or thickness), B. 
Pipes, et al. [10] defined the boundary-layer thickness as the distance 
from the edge, at which the interlaminar stress t is about 3 per cent of 
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the value calculated at the intersection of the ply interface and the edge of 
the laminate. The validity of this definition is somewhat questionable 
because the interlaminar stresses are singular at this point. In this 
study an alternative definition of the boundary-layer width is proposed on 
the basis of the strain energy density distribution in the composite laminate. 

The strain energy density distributions E(x,y) along the interface of 
the [6/-0/-0/0] graphite-epoxy laminates are shown in Fig. 8. The strain 
energy density remains relatively constant in the far field where classical 
lamination theory holds, and Increases drastically by an order of magnitude 
as the edge is approached. In this paper, the boundary- layer width B in a 
composite laminate is defined as the distance away from the edge where the 
strain energy density along the interface is three per cent higher than the 
nominal value obtained in the far field. In general, E(B,h'*’) differs 
slightly from E(B,h ) due to the discontinuous in-plane stress components 
at y = h"^ and h ; thus, an average value of B is designated as the width of 
the boundary-layer region. Based on this definition, values of B/W for 
[0/-0/-0/0] graphite-epoxy laminates are evaluated and shown in Fig. 9 
It is obvious from the figure that the [45°/-45°/-45°/45°] graphite-epoxy 
laminate has a higher value of B/W than those composites with other fiber 
orientations. As 0 changes towards either direction, B/W decreases rapidly. 
When 0 has an angle of [±0®] or [±90° ], B/W vanishes identically indicating 
that there is no boundary-layer effect in these cases since the two adjacent 
plies are identical. 

4.7 Effects of Fiber Orientation 

Effects of fiber orientation on the thermal boundary-layer response in 
[0/-0/-0/0] graphite-epoxy laminates are best illustrated by boundary-layer 
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thermal stresses along 0/-0 ply interface. The thermal stresses are observed 
to be significantly affected by the alternation of fiber orientation. Figures 
10 and 11 depict the distributions of thermal interlaminar shear stresses, t 

xy 

and y = h. The is relatively small for all fiber orientations 

studied, and the change of fiber orientation only alters the amplitude of the 

interlaminar shear stress slightly. The t reaches its maximum before it 

xy 

vanishes at the laminate boundary, where traction-free boundary conditions 
are satisfied exactly. The other thermal interlaminar shear component is 
more significant in the boundary-layer region than ^^^cy* ^yz ^ higher 

value along the interface of the [45®/-45°/-45®/45°] graphite-epoxy laminate 
than those in other ply ocnfigurations due to the higher value of 6^. As is 
expected, the interlaminar shear stress becomes unbounded as the laminate 
boundary is approached due to the stress singularity. As 0 moves towards 
either side, the amplitude of the interlaminar shear stress decreases. The 
distribution of interlaminar normal stress along the interface of angle-ply 
laminates with different 0 is shown in Fig. 12. The a is found to be 

y 

vanishingly small in the far field. It is small in tension first, then 
changes its sign, and becomes compressive as the free edge is approached. 

Like the interlaminar shear stress, becomes unbounded at the free edge (i-e. 
at X = b and y = h) and is significantly only within the boundary-layer region. 
In Fig. 13, distributions of in-plane thermal stress in [0/-0/-0/0] 
graphite-epoxy are given. The laminate thermoelasticity solution reveals 
that has a very small value in the far field, and remains relatively 
constant before rapid increase in its magnitude, as the laminate boundary 
is approached. This is apparently different from the prediction of classical 
lamination, theory, which suggests that, in the symmetric angle-ply 


[0/-0/-0/0] composite laminate, the in-plane thermal stress a vanishes 

z 
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throughout the composite. At the intersection of the free edge and ply 

interface, again has an unbounded value due to the stress singularity 

at that point. In the boundary-layer region, a along the interface in the 

z 

[±45]g laminate is found to have a higher value than in the composites with 
other ply configurations. As 0 changes from 45®, the magnitude of de- 
creases appreciably. In fact, the results shown in Fig. 13 Indicate that a 

z 

becomes negligibly small for the cases with 9 > 75°. 

4.8 Effects of Relative Ely Thickness 

Another important laminate variable investigated in this study is the 
effect of ply thickness or volume of the layer that is stressed inter- 
laminarly, since transverse deformation and failure in composites are affected 
significantly by lateral constraints in the laminates. It has been found 
that deformation and fracture in certain composite systems change significant- 
ly with the ply thickness while stacking sequence remains the same. In 
this paper, influences of ply thickness on the thermal boundary-layer stress 
in [45°/-45°/-45°/45°] graphite-epoxy laminate with various h^/W are examined 
(W being kept constant) . Following the aforementioned analytical procedure 
and solution scheme, a parametric study on thermal boundary-layer response in 
the composite laminates with various ply thicknesses k^/W has been conducted. 
Numerical results showing the effect of h^/W on thermal boundary-layer 
stress intensity factors are given in Table 3. The results indicate that, 
for a given laminate configuration and fiber orientation, the change of h^/W 
does not alter the thermal stress singularity but does affect thermal boundary- 
layer stress intensity factors appreciably. The composite laminate with 
h^/W - 0.4 0.5, seems to have smaller values of than other cases studied. 
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As h^/W approaches 0 or 1, the boundary-layer stress intensity factors reach 
their maxima. The higher values of developed in the composite by changing 
ply thickness h^/W provide an important basis for evaluation of initiation of 
transverse cracking and interply delamination in composite laminates. 
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5. SUMMARY AND CONCLUSIONS 

A study of thermal boundary-layer stresses in composite laminates has 
been presented. Formulation of the problem is based on the theory of 
anisotropic laminate thermoelasticity. With the aid of Lekhnitskii^s 
comp lex- variable stress functions, an eigenfunction expansion method is 
used to establish a system of coupled, governing partial differential 
equations for the problem. Numerical results for symmetric angle-ply 
[0/-e/-0/0] graphite-epoxy laminates are obtained. Effects of lamination 
and geometric variables on the thermal boundary-layer stress singularity 
and distributions are studied. Based on the information discussed in 
the previous sections, the following conclusions may be reached: 

1. Thermal stresses in the boundary-layer region of a composite 
laminate are inherently three dimensional in nature. They 
cannot be calculated by the classical lamination theory, but 
can be determined explicitly by the current approach. 

2. The thermal stress field in the boundary-layer region is 
singular in general. By using an eigenfunction expansion 
method, one can determine the order of the boundary-layer 
stress singularity by solving the characteristic transcendental 
equation obtained from the homogeneous solution. The order 

of the boundary-layer stress singularity depends on anisotropic 
thermoelastic properties of adjacent plies in the composite. 

3. The boundary-layer thermal stress field may be characterized by 
"thermal boundary-layer stress intensity factors" or "thermal 
free-edge stress intensity factors." The are functions of 
anisotropic thermoelastic constants of laminae, ply orientation 
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and laminate geometry, and may be used to evaluate strength 
degradation and initiation of interlaminar fracture 
(delamination) and transverse cracking. Their values can 
be determined by various methods such as the boundary 
collocation method and the finite element method. 

4. Thermal boundary-layer thickness which characterizes the domain 

where classical lamination theory does not hold can be 
determined explicitly by considering the change of strain 
energy density along the ply interface. The thermal boundary- 
layer thickness depends on lamination variables, geometric 
parameters, thermal loading conditions and thermoelastic ply 
properties. In graphite-epoxy composites, the case of 

0 = 45° possesses a higher thermal boundary-layer thickness 
than those of other fiber orientations studied. 

5. In comparison with previous approximate solutions, good 
agreement in the far field is observed, but appreciable 
discrepancy near the laminate boundary occurs. iThe difference 
is attributed to the fact that, in the boundary-layer region, 
the thermoelastic solution is completely governed by the 
singular terms which previous approximate solutions failed 

to include. 

6. The current method of approach is also valid for asymmetric 
composite laminates, since the bending, twisting and rotational 
components of deformation are included in the formulation. 
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Table 1 


First four nonzero eigenvalues for thermal boundary-layer stresses 
in symmetric angle-ply [0/-0/-0/0] graph! te-repoxy laminates* 


0° 




*2.3 


^4 

15° 

-0.64322 

E-3 

0.99670 

± 0.04191 

i 

1. 

30° 

-0.11658 

E-1 

0.95521 

± 0.15271 

i 

1. 

45° 

-0.25575 

E-1 

0.88147 

± 0.23401 

i 

1. 

60° 

-0.23346 

E-1 

0.83074 

± 0.27138 

i 

1. 

75° 

-0.89444 

E-2 

0.86469 

± 0.25007 

i 

1. 


corresponding to the strength of thermal boundary-layer stress 
singularity. 
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Table 2 

Thermal boundary-layer stress Intensity factors, K^, along the 
interface of [0/-0/-9/0] graphite-epoxy composite laminates*’*’ 

e K3 K5 Kg 

IS” 4.5996 E-1 -3.8534 E 2 -3.7120 E 1 -5.2624 E 3 8.0636 E 0 0 

30° 4.7031 E 0 -1.6466 E 2 8.6078 E 0 -5.7404 E 2 2.2357 El 0 

45° 1.4546 E 1 -1.9130 E 2 -7.4958 E-1 -4.1747 E 2 3.6662 El 0 

60° 1.3746 E 1 -1.7881 E 2 -2.2426 E 1 -3.7858 E 2 2.4394 El 0 

75° 4.5752 E 0 -1.5101 E 2 -2.9205 E 1 -4.9611 E 2 8.5717 E 0 0 

*hi = h2 = h, b = 8h 

’^’Values of are per °F change 
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Table 3 

Lamina thickness/volume effects on thermal boundary-layer stress 
intensity factors for [45°/-45'’/-45°/45®] Graphite-Epoxy Laminates*^ 



0.2 1,4870 E 1 -1.9554 E 2 -7.6620 E-1 -4,2672 E 2 3.7475 El 0 
0.4 1.4509 E 1 -1.9079 E 2 -7.4759 E-1 -4.1636 E 2 3.6565 El 0 
0.5 1.4548 E 1 -1.9130 E 2 -7.4958 E-1 -4.1747 E 2 3.6662 El 0 
0.6 1.4707 E 1 -1.9339 E 2 -7.5780 E-1 -4.2204 E 2 3.7064 El 0 
0.8 1.5580 E 1 -2.0487 E 2 -8.0275 E-1 -4.4708 E 2 3.9263 El 0 

*2b = 8W where W = constant = half laminate thickness 
"^Values of are per “F change 
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In-plane and Interlaminar Thermal Stresses along Interface in 
[45°/-45''/“45°/45''] Graphite-Epoxy Laminate, (h^=h 2 =h, b= 8 h) 


Through-Thickness Distribution of Boundary-Layer Thermal Stress 
in [45‘*/-45°/-45°/45° ] Graphite-Epoxy Composite (hjL=h 2 =h, b= 8 h) 


Through-Thickness Distribution of Boundary -Layer Thermal Stress 
Xy 2 in [45^/-45°/-45°/45*'] Graphite-Epoxy Composite 


Through-Thickness Distribution of Boundary-Layer Thermal Stress 
in [45°/-45®/-45°/45°] Graphite-Epoxy Composite 


Through -Thickness Distribution of Boundary -Layer Thermal Stress 
Oy in [45°/-45°/-45°/45°] Graphite-Epoxy Composite 


Strain Energy Density Distribution along Interface in [0/-0/-0/0] 
Graphite-Epoxy Composites (h^=h 2 =h, b/W= 6 ) 


Thermal Boundary-Layer Thickness versus Ply Orientation in 
[0/-0/-0/0 ] Graphite-Epoxy Composites 


Thermal Interlaminar Shear Stress t along Interface in [0/-0/-0/0] 
Graphite-Epoxy Composites (hjj^=fi 2 =h,^b= 8 h) 
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Graphite-Epoxy Composites ^ 
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cr^/AT (psi/^F) 


Fig. a Through-Thickness Distribution of Boundary-Layer Thermal 
Stress in [45“/-45V-45V45°] Graphite-Epoxy Composite 

(h^ =h2=h, b=8h) 
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Fig. 5 Through-Thickness Distribution of Boundary-Layer Thermal 
Stress in [45°/-45®/-45V45‘*] Graphite-Epoxy Composite 
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Fig. 6 Through-Thickness Distribution of Boundary-Layer 
Thermal Stress in [45‘»/-45®/-45«/45“] Graphite- 
Epoxy Composite 
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Fig. 7 Through-Thickness Distribution of Boundary-Layer 
Thermal Stress Oy in [45V-45‘’/-45®/45“3 Graphite- 
Epoxy Composite 
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Fig. II Thermal Interlaminar Shear Stress along Inter- 
face IN [0/-0/-0/0] Graphite-Epoxy Composites 
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Fig. 12 Thermal Interlaminar Normal Stress ay along Inter- 
face IN [e/-6/-0/e] Graphite-Epoxy Composites 
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Fig. 13 In-plane Thermal Stress along Int 
[ e/-0/-e/0] Graphite-Epoxy Composites 
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